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Abstract —In this paper, a receiver algorithm consisting of dif¬ 
ferential transmission and a two-stage detection for a single-input 
multiple-output (SIMO) phase-noise channels is studied. Specifi¬ 
cally, the phases of the QAM modulated data symbols are manip¬ 
ulated before transmission in order to make them more immune to 
the random rotational effects of phase noise. At the receiver, a two- 
stage detector is implemented, which first detects the amplitude 
of the transmitted symbols from a nonlinear combination of the 
received signal amplitudes. Then in the second stage, the detec¬ 
tor performs phase detection. The studied signaling method does 
not require transmission of any known symbols that act as pilots. 
Furthermore, no phase noise estimator (or a tracker) is needed 
at the receiver to compensate the effect of phase noise. This con¬ 
siderably reduces the complexity of the receiver structure. More¬ 
over, it is observed that the studied algorithm can be used for the 
setups where a common local oscillator or separate independent 
oscillators drive the radio-frequency circuitries connected to each 
antenna. Due to the differential encoding/decoding of the phase, 
weighted averaging can be employed at a multi-antenna receiver, 
allowing for phase noise suppression to leverage the large number 
of antennas. Hence, we observe that the performance improves 
by increasing the number of antennas, especially in the separate 
oscillator case. Further increasing the number of receive antennas 
results in a performance error floor, which is a function of the 
quality of the oscillator at the transmitter. 

Index Terms —Phase noise, multiple-antenna, phase averaging, 
distributed oscillators, Wiener process, differential modulation, 
two-stage detection. 

I. Introduction 

Phase noise due to phase and frequency instability in the local 
radio-frequency (RF) oscillators used in wireless communica¬ 
tion links results in synchronization issues, which degrade the 
system performance [1], [2]. These effects are more severe when 
high-order modulation schemes are used in order to attain high 
spectral efficiency [3]. It is also known that phase noise in RF 
oscillators increases with frequency [4], [5]. 

In multiple-antenna systems, the impact of phase noise is 
different depending on whether the RF circuitries connected to 
each antenna are driven by separate (independent) local oscilla¬ 
tors (SLO) or by a common local oscillator (CFO). Although the 
CFO configuration is of lower implementation complexity, the 
SLO configuration is unavoidable when a spacing, as large as a 
few meters, is needed between the antennas in order to exploit 


the available spatial degrees of freedom for higher multiplexing 
or diversity gains [6]-[10]. 

Carrier phase synchronization in multiple-input multiple- 
output (MIMO) systems has been extensively studied in the 
literature (e.g., [11]-[18] and the references therein). Phase syn¬ 
chronization methods in the case of the CFO configuration is 
similar to that in the case of a single-antenna system [l]-[3], 
which is followed by data detection that is generally used for 
MIMO systems [19]. However, designing receiver algorithms 
for joint phase noise estimation and data detection in the SLO 
case is more challenging. This is mainly due to the fact that 
first, coherent combining of the received signals at the receiver 
is not possible. Next, multiple phase noise parameters must 
be tracked and compensated before performing data detection. 
The effect of phase noise on the performance of MIMO sys¬ 
tems with a SLO configuration is studied in [9], [14], [18]. 
Authors in [15] investigate bounds on the mean-square error 
performance of phase noise estimators, including the extended 
Kalman filter (EKF). In [9], a joint phase noise estimation, data 
detection method based on maximum a posteriori (MAP) theory 
is proposed by applying the sum-product algorithm. However, 
most of the techniques proposed in prior studies suffer from 
high complexity, and are based on pilot-based transmissions for 
estimation and detection. 

In this work, we consider the scenario where a single-antenna 
user communicates on the uplink with a base station with mul¬ 
tiple antennas over an AWGN channel impaired by phase noise. 
For this scenario, we develop an algorithm that consists of dif¬ 
ferential transmission, followed by a two-stage amplitude and 
phase detector for detection of the transmitted signal (see [3] 
for two-stage detection). Differential transmission can be per¬ 
formed without the aid of pilot symbols, and no additional phase 
noise estimator/trackers are required at the receiver in order to 
track the random phase noise. Moreover, it has lower complexity 
than schemes, which are based on the MAP theory. More impor¬ 
tantly, due to the differential encoding/decoding of transmitted 
phase, weighted averaging can be employed at a multi-antenna 
receiver, allowing for phase noise suppression to leverage the 
large number of antennas. These aspects make the developed 
algorithm particularly attractive for systems i) that require low- 
complexity receiver algorithms, ii) whose performance are af- 


fected by the use of pilots, as in massive-antenna systems [20], 

[21] . We observe that the developed method is invariant to the 
oscillator configuration at the base station when applied to a 
SIMO phase noise channel. Our simulation results show that 
this method outperforms a Kalman filter-based algorithm [15] 
at high SNR or high phase noise scenarios. We also observe 
that when the SNR is high, the SLO yields a lower symbol 
error rate than the CLO because of the noise averaging effects. 
Furthermore, the performance of the studied method in SLO 
configuration improves by increasing the number of antennas 
till it achieves an error floor that is a function of the level 
of phase noise at the transmitter. This observation is in line 
with the results from prior work on the performance of uplink 
transmissions in presence of phase noise (see e.g., [10], [17], 

[ 22 ] ). 

Notations: Italic letters (x), boldface letters (x), and upper¬ 
case boldface letters (X) denote scalar variables, vectors, and 
matrices, respectively. The (i,i)th entry of matrix (X) is de¬ 
noted by [X]ij . With A/'(0, a^) and CA/'(0, cr^), we denote the 
probability distribution of a real Gaussian, and of a circularly 
symmetric complex Gaussian random variable with zero mean 
and variance cr^. For a given complex vector x, we denote by 
/X , the vector that contains the phase of the elements of x; 5ft( ), 
and 3(-), denote the real and imaginary parts of complex values, 
respectively; H-H denotes the Euclidean norm of a vector; the 
operator diag (•) generates a diagonal matrix from a vector; ly 
denotes the modified Bessel functions of the first kind with order 
v\ Q{-) is the Q-function. 


II. System Model 

Consider the uplink channel of a wireless communication 
system where a single-antenna user communicates with a base 
station equipped with M antennas over an AWGN channel 
impaired by phase noise. This yields the following 1 x M single¬ 
input multiple-output (SIMO) phase-noise channel: 

yfe = ©fehxfe-I-Wfe, fc = l,...,n, (1) 


where Xk is the transmitted symbol in the kth time instant, 
0fc = ,... and em,k = + 6^^, for 

m = 1,..., M, where and represent the random 

phase noise samples, in the fcth time instant, from the local os¬ 
cillators at the transmitter and receiver, respectively. We model 
the phase noise as a Wiener process [23] 


,[t] 


fc 



= 

= el 


[t] 

fc-i 


,fc-i 


jt] 

^fc-i 

-a" 


fc-l! 


( 2 ) 

(3) 


Wfc = ... ,WM,kY' are the AWGN samples, which are 

drawn independently from a CAf(0,2) distribution.^ 

III. SLO Configuration 

In the SLO configuration, coherent combining of the received 
signals is not possible. In this section, we study a differential 
transmission algorithm, which is followed by a two-stage data 
detection scheme at the receiver. The studied signaling method 
can be used along with the various modulation schemes such as 
M-QAM. This method enables a non-data-aided carrier phase 
synchronization at the receiver over the SIMO phase-noise chan¬ 
nel. 

As illustrated in Fig. 1, at the transmitter side, we first mod¬ 
ulate the data bits as QAM symbols denoted as Sfc, which in 
the polar coordinate can be expressed as Sk = Vke^'^^. Next, 
we manipulate the phase of each symbol before sending it over 
the SIMO phase noise channel. More specifically, the phase of 
the transmitted symbol is generated by accumulating the phase 
of the current symbol and that of the previously transmitted 
symbol, i.e., 

Xfe = (4) 

Note that the amplitudes of the symbols remain unaltered. 

At the receiver side, we propose a two-stage detection proce¬ 
dure in order to detect the transmitted QAM symbols. First, a 
symbol-by-symbol maximum a posteriori (MAP) detection of 
the amplitude is performed. Then the transmitted signal ampli¬ 
tude, as determined by the amplitude detector, is used along with 
the phase of the received signals at each antenna in order to 
perform differential maximum likelihood (ML) phase detection. 
It is straightforward to show that given Tfc, tk = ||yfc|p follows a 
noncentral chi-squared distribution with 2M degrees of freedom 
and noncentrality parameter r^,||h|p [10], i.e., 

M 

tk ^ \rk\\h\\ + Wi\^ +'^ \wm\^ , (5) 

m—2 

where ~ denotes equality in distribution, and {wm}, m = 
1,... ,M are independently and identically distributed (i.i.d.) 
from CAf{0, 2). The MAP amplitude detector is determined as 

ffe = argmax f{rk \ tk) = argmax f{tk \ Tfc)p(rfc), (6) 

rk rk 

where p{rk) denotes the probability mass function of rk- The 
likelihood function f{tk \ rk) can be determined from the defi¬ 
nition of a noncentral chi-squared distribution [26, Eq. 2.44] 


where and {A|f^} are the Gaussian random samples 

drawn independently from Af{0, ) and Af{0, o'? ) distri¬ 

butions, respectively'. The vector h = [hi,..., hMf contains 
the path-loss coefficients, which are assumed to be deterministic, 
time-invariant, and known to the receiver. Einally, the entries of 


*For discussions on the limitations of this model see [2], [24], [25] and the 
references therein. 


f{tk\rk) = 

1 / i \ (AT-l)/2 

y ^M-t(r.||h||V4).(7) 

The detector can be implemented by substituting (7) in (6), and 
then finding the maximizer by evaluating (6) over all possible 

^As we shall see, normalizing the noise variance to 2 will turn out convenient. 




Fig. 1. System diagram of the proposed differential signaling method for an 1 x M SIMO channel. 


symbol amplitudes. For example, there are three amplitude pos¬ 
sibilities when using 16-QAM, 
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on the phase of the transmitted symbols. From (13) we observe 
that /ym,k — /Vm,k-i, for m = are independent, 

and approximately follow a Gaussian distribution. Hence, the 
sufficient statistic for detection of <f>k is given by [27] 

1 ^ 

= i /ym,k - /yra,k-l ) ■ (14) 

m—1 

It is straightforward to show from (13) and (14) that 


where E denotes the average symbol energy. 

In the next stage, a differential phase detection is performed 
as follows. The phase of the received signal per each antenna is 
written as [14, Eq. 12] 
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where Wm,k = ^{wm,k}/{fk \hm\) is a zero-mean real Gaus¬ 
sian random variable. The approximation in (11) is valid at 
moderate and high SNR [14]. In the first step of phase detection, 
from the phase of the current received signal, we subtract the 
phase of the previously received signal at each receive antenna 
as 


As we see in (16), in the SLO case, the total variance of phase 
noise from receiver’s local oscillators reduces upon increasing 
M. This is because the receiver phase-noise samples are i.i.d. 
over receive antennas. ^ 

In order to detect (f>k we perform a maximum likelihood 
detection as follows; 


/ym,k /ym,k — l ^ /Xk /Xk — 1 ~t~ Om,k ^m,k—l 

4“ '^m,k '^m,k — l (12) 

= 4:k + aW . + aW . 

4” '^m,k '^ra.k — l: (15) 

where the equality in the last step is obtained by using (11), (4), 
(2), and (3). By means of this differential transformation, we 
mitigate the effect of instantaneous unknown phase noise values 


(j)k= argmax /(V'fe I (/ife) (17) 

^ke /C{fk) 

= argmin {'ipk - (t>kf, (18) 

<l>k& /C{fk) 

^Note that (14) is an equally-weighted sum of phases over the receive anten¬ 
nas. However, the weights can be optimized in order to deal with differences 
in quality of the signals received at the different antennas. Similaiiy, the ampli¬ 
tude detection can be improved by considering an optimized weighted sum of 
amplitude squares of received signals. 









































































where the minimization in (18) is performed over the phase of 
a set of transmitted QAM symbols with amplitude fk denoted 
as /C{fk). Note that the phase detection in (18) can be seen as 
a Euclidean distance detector in the phase domain. 

IV. CEO Configuration 

In the CEO configuration, we have that 6i^k = • • • = 0m, fc = 
dk for all k. Therefore, the input-output relation (1) simplifies 
to 

yk = e^^>‘hxk + Wk. (19) 


depends on amplitude of the current symbol and that of the pre¬ 
viously transmitted symbol. In order to simplify the analysis and 
remove the dependency on the past symbol, we compute cr^j, by 
setting rk-i = E (the average symbol energy). 

Error floor at liigh-SNR, when M oo: In order to 
determine the SEP at high-SNR for large number of antennas, 
the pair-wise probability of error can be simplified by evaluating 

lim = q( , (23) 

M-s-oo / 


As the M signals received in a given time instant are affected 
by only one phase noise process, it is possible to first coherently 
combine the received signals by projecting on h/||h||, and 
convert (19) to a single-input single-output (SISO) channel, fol¬ 
lowed by the compensation of Ok and data detection. Note that 
the CEO case, after coherent combining, is a special case of the 
SLO configuration with M = 1. 

V. Symbol Error Probability Analysis 

In this section, we derive analytical expressions for the sym¬ 
bol error probability (SEP) performance of the developed dif¬ 
ferential algorithm and the two-stage detector for arbitrary con¬ 
stellations. Consider the transmission of an arbitrary constella¬ 
tion with N equally-likely symbol points in (1). The symbol 
error probability for this constellation is upper-bounded by 
averaging over all pair-wise symbol error probabilities (union 
bound) [28] as 

^ N N 

E (20) 

where P{Eij) is the pair-wise symbol error probability. This 
corresponds to the probability of the error event Eij, where the 
symbol j is detected given that symbol i has been transmitted. 
Eor a two-stage detector, the pair-wise symbol error is a result 
of two events: the error in pair-wise detection of the amplitudes, 
denoted as and the error in detection of the phase of 

^ phase 

the transmitted symbol, denoted as E^^ . Accordingly, the pair¬ 

wise symbol error probability becomes 

P{E,j) = p(£;™Pl“^‘le^p^^phase^_ ^21) 

In order to simplify the analysis, we consider a high- 
SNR scenario, where no amplitude detection error occurs, 
i.e., = 0. In this case, the error in sym¬ 

bol detection occurs only in the phase domain. Therefore, 

= 0 for symbols i and j with different ampli¬ 
tudes, while p(p™P*““‘^‘=) = X for equal-amplitude symbols. 
By using (15), the pair-wise phase detection error probability 
for equal-amplitude symbols can be found as 

P{E^^^) = Q , (22) 

where (f)\, and (jy\, denote the phase of symbols i and j, respec¬ 
tively. Note that the value of , which is computed from (16) 
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We can observe from (24) that the error floor in the CEO case is 
a function of phase noise innovation variance in transmitter and 
receiver oscillators, while in the case of SLO configuration, it is 
only a function of phase noise innovation variance of transmitter 
oscillator. 


VI. Overview oe the Kalman Eilter Method 


In this section we present a EKE-based phase noise 
estimation-data detection algorithm similar to that of developed 
in [15]. The QAM symbols are transmitted directly over the 
phase-noise channel (1) while an EKE at the receiver compen¬ 
sates the effect of phase noise before data detection. Note that 
the QAM symbols transmitted here have not been manipulated 
as in (4). 

As already mentioned, in the SLO configuration, coherent 
combining of the received signals is not possible, and the vector 
Sk = ; 0M,fc]^ consisting of M phase noise processes 

must be tracked. The state equation in this case is 


— 0fc-i + ^k-ii (25) 


where Aj, ^ A/^(0, S) and 
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for i, j G {1,..., M}. We use (26) along with the observation 
model (1) to estimate 6k by using an EKE [27]. Note that when 
pilot symbols are not transmitted, the EKE does not have reliable 
estimates of Xk, thereby causing the performance of the EKE 
to deteriorate. In order to avoid this problem, we transmit pilot 
symbols, which are known at the receiver. More specifically, at 
each time instant, we first perform the prediction step of the EKE 
which is independent of the observation. If a pilot symbol is 
transmitted, we use the pilot transmission to perform the update 
step. Eor the data symbols, we first use the predicted phases to 
de-rotate the received signal at each antenna. Then a Euclidean- 
distance based detector is used to detect the transmitted sym¬ 
bol [14]. Einally, the detected symbol is used in the update step 
of the EKE. 






Fig. 2. SEP versus SNR for a 1X 10 SIMO channel with the CEO configuration. 



Eig. 3. SEP versus SNR for CEO and SLO configurations. In this simulation 
= 0.01 [rad^] and M = 10. 

VII. Results and Discussions 

In this section, we present the performance results of the pro¬ 
posed signaling method using 16-QAM in terms of SEP versus 
SNR per symbol (E/2). We scale down the transmit power by 
M when increasing the number of receive antennas to keep the 
average receive SNR constant. In our simulations we consider 
a quasi-static fading model [29, p. 2631], [30, Sec. 5.4.1] for 
the channel, where the {hm} are independently drawn from 
aCM{0, 1) distribution. The channel is kept a constant over each 
simulation trial, and is known to the receiver. The variance of the 
phase noise innovation for transmitter and receiver oscillators 
are set equal, unless it is explicitly mentioned. 

Fig. 2 compares the SEP of the proposed differential method 
(DIF), versus the EKE algorithm for a 1 x 10 SIMO phase-noise 
channel with CEO configuration for various SNR and cr^ values. 



Fig. 4. Analytical and simulated SEP versus number of antennas and various 
phase noise innovation variances. In this simulation SNR = 40 IdB]. 


For the EKE algorithm, a pilot symbol is transmitted every 50 
data symbols resulting in a pilot density of 2%, while for the 
differential method no pilot symbols are transmitted. When the 
phase noise innovation variance is high, e.g., = 0.01 [rad^] 

or = 0.005 [rad^], the differential method performs sig¬ 
nificantly better than the EKE algorithm. The EKE method 
performs close to the differential method for a\ = 0.001 for 
moderate and high values of SNR. However, as the EKE is 
highly sensitive to the errors in data detection (i.e., when the 
transmitted symbol is not a pilot symbol), at low SNR, the 
performance of the EKE is much worse than the differential 
method even for low values of a\. 

Fig. 3 illustrates the SEP versus SNR of the proposed algo¬ 
rithms for the CEO and the SLO configurations when M = 10. 
It can be seen that the EKE performs significantly better in the 
SLO configuration compared to the CEO setup. Similar results 
can be observed for the differential method at high SNR. This 
gain can be explained as follows: in the SLO case, we receive M 
independent noisy observations of the phase of the transmitted 
signal, which enables us to perform an averaging before detec¬ 
tion. As can be seen in (16), this manifests as a reduction of 
the effective variance of the noise that distorts the phase of the 
signal. Similar observations have been reported in the capacity 
calculations for the SIMO phase noise channel [10]. However 
at low SNR, the proposed differential method performs worse 
in the SLO configuration. This is due to the fact that SNR per 
antenna becomes quite low as it is scaled by a factor of M = 10, 
which renders the approximation in (11) to be erroneous. The 
high-SNR approximation of the SEP performance of the differ¬ 
ential signaling method is also presented in Fig. 3. It can be seen 
that the analytical SEP matches the simulation results for SNR 
values around 30 [dB] and above. 

Fig. 4 compares the analytical and the simulated SEP of the 
SIMO system with SLO configuration, for various values of 
phase noise innovation variance at the transmitter, versus the 




























































number of antennas. An error floor is observed by increasing M. 
We see that the error floor is a function of the phase noise at the 
transmitter. 

VIII. CONCLUSIONS 

In this paper, we developed a receiver algorithm for a SIMO 
phase-noise channel that consists of a differential signaling 
method combined with a two-stage detection of the transmitted 
signals. This method can be used in strong phase noise scenarios, 
and when the pilot density mandated by the application is low. 
It can also be useful in massive-antenna systems because of its 
low complexity compared to conventional phase noise tracking 
algorithms like the EKF, and also in the systems whose perfor¬ 
mance is limited by pilot contamination. The performance of 
the developed receiver improves by increasing the number of 
receive antennas due to noise averaging effects. Specially when 
independent oscillators are used at the receiver, by increasing 
the number of receive antennas, the effect of phase noise from 
receiver oscillators averages out. In this configuration, further 
increasing the number of antennas results in an error floor in the 
performance, which is a function of the quality of the oscillator 
used at the transmitter. 
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